Coherent backscattering of Bose-Einstein condensates in two-dimensional disorder 

potentials 
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We study quantum transport of an interacting Bose-Einstein condensate in a two-dimensional 
disorder potential. In the limit of vanishing atom-atom interaction, a sharp cone in the angle- 
resolved density of the scattered matter wave is observed, arising from constructive interference 
between amplitudes propagating along reversed scattering paths. Weak interaction transforms this 
coherent backscattering peak into a pronounced dip, indicating destructive instead of constructive 
interference. We reproduce this result, obtained from the numerical integration of the Gross-Pita- 
evskii equation, by a diagrammatic theory of weak localization in presence of a nonlinearity. 
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The past years have witnessed an increasing number 
of theoretical and experimental research activities on the 
behaviour of ultracold atoms in magnetic or optical d.\s- 

order potentials m [a [31 n E 1 E la 1 uni nn ini [ni im • 

A central aim in this context is the realization and un- 
ambiguous identification of strong Anderson localization 
with Bose-Einstein condensates, which was attempted by 
several experimental groups [D |5J |3] with recent success 
[H [S] , and theoretically studied both from the perspec- 
tive of the expansion process of the condensate [6^ as well 
as from the scattering perspective [71 [5]. Complementary 
studies were focused on localization properties of Bogoli- 
ubov quasiparticles [HI [12] , on dipole oscillations in pres- 
ence of disorder (TT] [T^j , as well as on the realization of 
Bose glass phases [I31IT4]. 

The above-mentioned topics mainly refer to processes 
that are essentially one-dimensional (ID) by nature. 
Qualitatively new phenomena, however, do arise in two 
or three spatial dimensions, due to the scenario of weak 
localization. The latter manifests in a slight reduction of 
the transmission probability of an incident wave through 
a disordered region as compared to the classically ex- 
pected value, due to constructive interference between 
backscattered paths and their time-reversed counter- 
parts. This interference phenomenon particularly leads 
to a cone-shaped enhancement of the backscattering cur- 
rent in the direction reverse to the incident beam, which 
was indeed observed [TS] and theoretically analyzed [T5] 
in light scattering processes from disordered media. Re- 
lated weak localization effects also arise in electronic 
mesoscopic physics, leading to characteristic peaks in the 
magneto-resistance [17l[T8]. 

In this Letter, we investigate the phenomenon of co- 
herent backscattering with atomic Bose-Einstein conden- 
sates that propagate in presence of two-dimensional (2D) 
disorder potentials. An essential ingredient that comes 
into play here is the interaction between the atoms of the 
condensate. On the mean-field level, this is accounted for 
by the nonlinear term in the Gross-Pitaevskii equation 



describing the time evolution of the condensate wavefunc- 
tion. Indeed, nonlinearities do also appear in scattering 
processes of light e.g. from a gas of cold atoms, due to 
the saturation of the intra-atomic transition [191 HOI [21] • 
In this case, however, the saturation also leads to inelas- 
tic scattering [SU] [2T] and, in addition, the nonlinear- 
ity competes with other dephasing mechanisms induced, 
e.g., by polarization phenomena [55] or thermal motion 
[55] . The complementary process of atomic condensates 
scattering from optical random potentials in the mean- 
field regime provides a cleaner situation where the coher- 
ence of the atomic wavefunction remains well preserved 
in the presence of the nonlinearity. As we shall argue be- 
low, this leads to substantial modifications of the coher- 
ent backscattering feature. In particular, the interaction 
turns constructive into destructive interference, leading 
to a negative coherent backscattering peak height. This 
is reminiscent of the weak antilocalization effects due to 
spin-orbit interaction observed for mesoscopic magneto- 
transport [24] . 

The starting point of our investigation is the time- 
dependent 2D Gross-Pitaevskii equation describing the 
mean-field dynamics of the condensate in presence of the 
disorder potential V{r) [r = (a;,j/)], 

ih^^i^{r,t) = (-!^^ + V{r) + ~g{x)\i>{r,t)Ai:{f,t) 



S{t)5{x — Xq) exp{—i^t/h) , 



(1) 



where S(t) denotes a source term simulating the coherent 
injection of matter waves with chemical potential /i from 
an external reservoir onto the scattering region [7]. In 
the numerical integration of Eq. (IT]), S{t) is adiabatically 
increased from zero to a final value 5*0 that corresponds 
to a fixed incident current density jin- Periodic boundary 
conditions are imposed on the transverse boundaries (in y 
direction) of the numerical grid to ensure a homogeneous 
flow in absence of disorder, whereas absorbing boundary 
conditions applied at the edges of the longitudinal {x) 
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FIG. 1: (Color online) Scattering geometry and stationary 
scattering state associated with a randomly generated disor- 
der potential. The left-hand side displays V{x,y) in a gray- 
scale plot and shows the spatial variation of the nonlinear- 
ity g{x). The upper right panel shows the density of the 
corresponding scattering state that is populated through the 
numerical integration of the inhomogeneous Gross-Pitaevskii 
equation |T]). The lower right panel shows the decay of the 
coherent mode KV")!^ and the density (IV'P) with x, averaged 
over y for ~ 10^ randomly generated disorder configurations. 
Parameters: kL = 40, kW = 120, kcr = 0.5, Vo = 0.614p, 
g = 0.005, Jin ~ hk^ /m, with k = y^2mii/h. 



direction allow us to inhibit artificial backreflection of 
outgoing waves with rather high accuracy [25j . 

In Eq. (Ilj , the effective 2D interaction strength is writ- 
ten as g{x) = h? g{x) / {2m) , with the dimensionless non- 
Hnearity parameter g{x). In presence of a harmonic con- 
finement of the condensate in the third spatial dimen- 
sion with the oscillator length a±{x) = ^Jh/[mbj±{xy\, 
we have g{x) = A^/2^Tas / a±{x) , where a^ denotes the s- 
wave scattering length of the atoms. We assume that 
g{x) is adiabatically ramped on and off in front of and 
behind the disorder region, as shown in Fig. [l] Physi- 
cally, this spatial variation of the nonlinearity, which is 
needed in order to avoid nonlinear effects at the posi- 
tion of the source and the absorbing boundaries, would 
correspond to a finite extent of the transverse harmonic 
confinement into which the condensate is propagating. 
As for the disorder potential V{r), we choose a Gaus- 
sian random process characterized by a vanishing mean 
value {V{r)) = and a Gaussian correlation function 
{V{f)V{r + Ar)) = Vg^e^'^'''/^'^' with correlation length 
a. We focus in the following on the parameters ka — 0.5, 
with k = ^j2/m^/h the wavcnumber of the incident beam, 
and Vq/ ^ — 0.614. The incident current density reads 
Jin = ^fc|'0oP/'7i, where we set ipQ = k for the amplitude 
of the incident wave p^ . 

At the above values for a and Vq, scattering in the 
disorder region is approximately isotropic. This is quan- 
titatively expressed by the equivalence of the two rele- 
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FIG. 2: Scattering mean free path Is (left panel) and trans- 
port mean free path Itr (right panel) in the disorder poten- 
tial for ka = 0.5 in absence of the nonlinearity. In the left 
panel, the numerically computed lengths are compared with 
the Born approximation S (solid line) and in the right panel 
with the Boltzmann mean free path ([3| (^s, dashed line) and 
the expression Q that takes into account weak localization 
corrections {Itr, solid line). We find Itr — is for Vo < p, 
which characterizes isotropic scattering. 



vant length scales that the disorder averages introduces 
for the transport process of the condensate: the scatter- 
ing mean free path £s , which describes the average decay 
of the incident coherent mode inside the disorder region 
according to \{ip{r))\'^ oc exp{—x/£s), and the transport 
mean free path itr, which characterizes the decay of the 
average density (|V'(^)P) (see Fig. [l]). In absence of the 
nonlinearity, the scattering mean free path is in leading 
order in Vq given by the Born approximation 



{kQ- 



{7T/2){Vo/fi)'{kaflo{k^a^) cxp{-k'a') (2) 



where lj{0 is the modified Bessel function of order j. 

The transport mean free path can be extracted from 
the linear decrease of (|'0(r)P) with x according to 
(|V'(?OP) (X {L + zoitr ~ x), with zo = 0.82 in two spatial 
dimensions, and L the longitudinal extent of the disorder 
region. In lowest order in Vq, £tr is given by the Boltz- 
mann transport mean free path £b defined through 



',2_2\ 



ejeB^l-h{k'cT')/Io{k'a') 



(3) 



Weak localization effects lead to logarithmic corrections 
that yield for kls > U2Z1 HI] 



£B[l-2(Ms)-4og(L/£B)]. 



(4) 



As shown in Fig. [2] the expressions Q and (El are in 
good agreement with the numerically computed values 
of is and itr for Vq < fi. Specifically at ka = 0.5 and 
Vo/n = 0.614, we find fc4 ^ 9.61 and kitr ^ 9.75. 

The angle-resolved current in backward direction is nu- 
merically computed from the decomposition of the re- 
flected wave 'ip,-ci{x,y) = 'ip{x,y) — ipQexp{ikx) at fixed 
position X close to xq [where g{x) is negligibly small] into 
the transverse eigenmodes Xn{y) '^ exp{imTy/W), which 
support outgoing waves into the directions with the an- 
gles 9n = arcsin[27rn/(A:VF)]. Figure Is] shows the average 
angular density j{9) of the backscattered current, which 
is normalized such that L j{9)dd = 27r. In the linear 
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FIG. 3: (Color online) Angle-resolved current density of 
backscattered atoms in absence and in presence of the non- 
linearity, obtained from the average over ~ 10^ disorder con- 
figurations (parameters as in Fig. [I] the error bars denote the 
statistical standard deviation). The coherent backscattering 
cone for g = (black line) is transformed into a pronounced 
dip for intermediate nonlinearities {g = 0.02, bold red line), 
and turns into a smooth peak structure at larger values of 
g {g = 0.06, dashed blue line). The inset shows the angle- 
resolved current for the case of a tilted incident beam where 
the source term in Eq. ([l| populates the transverse eigen- 
mode defined by the angle 66 — 0.32. In contrast to the 
smooth peak, the cone and dip structures are indeed found at 
the angle that corresponds to retro-refiection of the incident 
beam, which confirms that they both arise due to interference 
between refiected paths. 



case {g = 0), we encounter the well-known cone struc- 
ture at 6* = 0, which is a characteristic signature of weak 
localization [T21[TB]. Rather small values of g ~ 0.02 cor- 
responding to g\ip{r)\'^ ~ 10~^/J,, are sufficient to substan- 
tially modify this cone-shaped peak. Most interestingly, 
it is not washed out by the nonlinearity, but transformed 
into a dip that roughly has the same shape as the peak 
at (7 = 0. This indicates that the underlying interference 
phenomenon between reflected scattering paths is still 
effective at finite g, but has turned from constructive to 
destructive. 

The occurence of a dip in the backscattered current is 
confirmed by calculations based on the diagrammatic ap- 
proach for weak localization in presence of a nonlinearity 
pn [50] . Assuming the realization of a stationary scat- 
tering state, the average density (|V'(^P) is expressed 
in terms of ladder diagrams, which amounts to neglect- 
ing interference, and thus describing wave transport as 
a classical random walk. This assumption is valid ap- 
proximately for a dilute medium, i.e. for k£ ^ 1 with 
£ = £b — ^tr — ^s (for isotropic scattering). Further- 
more, we assume the condition g'^\ipo/k\'^k£ <^ 1 under 
which scattering from the fluctuations 5|'0(OP of the 
nonlinear refractive index is negligible compared to scat- 
tering from the disorder potential V{f) [29]. Therefore, 
the average density (\'4'{^\'^) remains approximately un- 
affected by the nonlinearity, and thus is well described 
by linear transport theory. From the average density. 
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FIG. 4: Backscattered current at 6 = as a function of the 
nonlinearity g (parameters as in Fig. [l]), obtained from the 
numerical simulation (symbols) and from the diagrammatic 
theory, Eqs. ([Mt]) (solid line). The horizontal dashed line 
indicates the diffuse background intensity jiL(O). Negative 
cone heights jc(0) < leading to a dip in the angle-resolved 
current density j{d) appear for g > 0.01. 



the flux backscattered in direction 6 = results as 

In a second step, the coherent backscattering peak 
is calculated by means of crossed (Cooperon) dia- 
grams, describing interference between reversed scatter- 
ing paths. Following the diagrammatic approach pre- 
sented in Ref. [SU] , we obtain the height of the coherent 
backscattering peak from the transport equations 



Ci{x) 



\^^o\^e-'^'[l + -J dx'gix')C^{x')], (5) 



^ dx' 

^ 



Kn 



X ~ X 



(Ci(x') + C,(a;')) + 



d^MI' 



X 



(6) 



r'niax(a;,a:') 

X / dx"g{x")(Ci{x") + C,{x")) 

inin(x,a;') 



for the "Cooperon intensity" Ci{x) and the "coherent 
Cooperon intensity" Cc{x), with x = max(a;,0) and Kqi 
the modiflcd Bessel functions of the second kind. The 
contribution to the flux scattered in backward direction 
then results as 



"■^ dx 



-(l^(x)l^ 



''/'(c,ix)+ 
dx'g{x')Ci{x') 



(7) 



Note that nonlinear processes also occur for xq < x < 
where V{r) = but g{x) > (see Fig. [l]). Hence, the 
cone height jc(0) — in contrast to the background inten- 
sity jl(0) — explicitly depends on the spatial extent of 
the nonlinearity region in front of the disorder potential, 
and can therefore be tuned through the ramp-up of g{x). 
In the absence of nonlinearity, the above equations re- 
duce to linear transport theory in a two-dimensional slab, 
yielding C,{x)\^^, = (IV'(a^)P) - |VoPexp(-a:/£). We 



then obtain jc(0)L^o = Jl(0) — 1/2, which expresses reci- 
procity symmetry, i.e. the equality of reversed path am- 
plitudes (the term 1/2 describes single scattering). For 
g j^ 0, however, the nonlinearity turns Ci{x) into a com- 
plex quantity, as evident from the terms proportional to 
ig in Eqs. (5]|7l. This indicates an effective phase dif- 
ference between the reversed scattering paths. Conse- 
quently, the backscattered current jc(0) is expected to 
decrease with increasing nonlinearity, and may even be- 
come negative if this phase difference is sufficiently large. 

This latter situation is indeed encountered if the set 
of equations (IsIIt]) is numerically solved for the system 
parameters under consideration. As shown in Fig. [4] the 
total flux jl(0) -|- jc(0) resulting from Eqs. (5][7l (solid 
line) agrees rather well with the average value for j{0) 
obtained from the numerical simulation (symbols). Dis- 
crepancies are attributed to weak localization corrections 
in the background intensity, which would specifically lead 
to a reduction of the backscattered flux a,t g — 0, and 
to an additional contribution to the Cooperon intensity, 
termed C2{x) in Rcf. [50], which was neglected in the 
derivation of the above transport equations. Details on 
these additional ingredients will be presented elsewhere. 

At larger nonlinearities, g > 0.03, the numerical prop- 
agation of the inhomogeneous Gross-Pitaevskii equation 
(fTl does not converge to a stationary scattering state, 
but leads to a permanently time-dependent behaviour of 
ip{r,t), as predicted in Refs. [53] and encountered also 
in the transport of condensates through ID disorder po- 
tentials [3. In this regime, the average backscattered 
current again displays a peak around 9 — 0; this peak is, 
however, comparatively broad and does not arise from 
a coherent backscattering phenomenon. This becomes 
obvious if we inject the incident wave with a finite angle 
9q ~ 0.32 (corresponding to the transverse eigenmode xe) 
onto the disorder region. While the cone and dip struc- 
tures a,t g = and 0.02 appear, as shown in the inset of 
Fig. [3] at the expected angle of coherent backscattering, 
corresponding to retro-reflection of the incident beam, 
the broad peak at g = 0.06 is not affected in this way. 

In conclusion, the presence of a small nonlinearity re- 
verts the scenario of weak localization and gives rise to a 
cone-shaped dip, instead of a peak, in the angle-resolved 
backscattered current density. This phenomenon appears 
to be rather robust; it is numerically encountered also for 
disorder potentials with longer correlation lengths a giv- 
ing rise to anisotropic scattering, and we expect its man- 
ifestation also in three spatial dimensions (as predicted 
by the diagrammatic theory) as well as for speckle dis- 
order where diagrammatic approaches would have to be 
based on the treatment of Ref. [3S]. We therefore believe 
that the effect would be measurable, for a reasonably 
large range of parameters, in state-of-the-art transport 
experiments with coherent Bose-Einstein condensates in 
well-controlled disorder potentials. 
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